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Abstract 

We examine the quantum energy levels of rectangular billiards with a point- 
like scatterer in one and two dimensions. By varying the location and the 
strength of the scatterer, we systematically find diabolical degeneracies among 
various levels. The associated Berry phase is illustrated, and the existence of 
localized wave functions is pointed out. In one dimension, even the ground 
state is shown to display the sign reversal with a mechanism to circumvent 

the Sturm-Liouville theorem. 
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A sense of mystery has accompanied the geometric phase from the day of its discovery 
It has deepened when the connection between the geometric phase and the gauge 
field has been recognized ||. Although there are several model systems in relatively simple 
settings, few qualify as transparent examples that give us the intuition into the origin and 
the workings of the geometric phase. 

In this Letter, we examine the candidates for just such examples. We study the bound 
state problem in a rectangular billiard with a pointlike scatterer inside. The zero-size scat- 
terer is an idealized limit of a small obstacle. In this limit, the solution of the problem is 
expected to be simplified substantially. That simplification is indeed realized when the di- 
vergence caused by the singularity of the interaction is handled properly with the self-adjoint 
extension theory of functional analysis 0-0. We look at the energy level structure of this 
system along with its one-dimensional analogue. We point out the existence of a set of dia- 
bolical points and of associated geometric phase. An examination of wave functions around 
the diabolical points also reveals the emergence of unexpectedly localized wave functions. 

We first consider a one-dimensional delta potential problem with the Dirichlet boundary 
condition (vanishing wave function) at x = and 1. We place a pointlike scatterer of 
strength 1/s at x = x : 

V = - 8(x-x ). (1) 
s 

The diagonal element of the free Green's function at x = Xq is given by 

GH = £ (2) 

n=l U ~ £ n 

where e n = n/A ■ n 2 and (p n (x) = V% sin(ri7rx) (n = 1,2, ...) are the eigenvalues and eigen- 
functions of the unperturbed system. (We chose the mass to be 2tc.) The eigenvalues of the 
problem with the delta function are given as the solutions of the equation 

G(u) -8 = 0. (3) 

The eigenfunction (apart from the normalization) corresponding to an eigenvalue u a is given 
by 



oo 



_2 CUq £ n . 




(4) 



In Fig.l, we plot the first six eigenvalues at several function of s. Fig. 1 fa) 

shows the existence of level crossings at infinite strength s = with xq = 0.5. Fig. 1 (b) 
displays avoided crossings around s = for xq = 0.47, indicating that the crossings at Xo = 
0.5 are diabolical. That this is indeed the case is proved by examining Fig. 2 (a), where we 
display wave functions at various values of (xq, s) around the point (0.5, 0) starting from the 
ground state at (0.5, 0.1). By going round the point (0.5, 0), the wave function experiences 
the Herzberg-Longuet-Higgins sign reversal [H], which is a special case (co-dimension two) 
of the Berry's geometric phase for real symmetric Hamiltonians. The mechanism of the sign 
reversal in this simplest of all examples is easy to comprehend: It is a consequence of two 
conflicting ways of connecting nodeless and one-node wave functions. 

At first sight, this result seems to contradict with the Sturm-Liouville theorem of no- 
crossing and no- node |§, which states that in one dimension there is no level crossing and 
the ground state always stays nodeless. However, at the diabolical point (xq,s) = (0.5,0), 
the strength of the potential is infinite, where the no-crossing theorem does not have to hold. 
Also, as seen in Fig. 1, the ground state at s > side is smoothly connected to the second 
state with one node at s < side. This happens because when one goes from positive to 
negative s, a disconnected level appears from the negative infinity to become the ground 
state for s < 0. In this rather tricky manner, the ground state of s > side is permitted to 
acquire the sign reversal through the rotation around the diabolical point. 

Looking at Fig. 1 (a), we notice that there are a series of pairwise degeneracies at 
(xq,s) = (0.5,0). This occurs between the states with even n which are unaffected by the 
delta potential at the center and the states with odd n which lie next to them. Similarly, one 
realizes from Fig. 1 (c)-(e) that the third state, sixth state, ninth state, etc. get crossed by 
the states next to them at (xq, s) = (1/3, 0). More generally, if one assumes that N and M 
are the relative primes satisfying N > M > 0, the location (xq, s) = (M/N, 0) is a diabolical 
point for the (N — l)-th and N-th states, (2N — l)-th and 2iV-th states, etc. (all counted 
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at s > side). Clearly, one exhausts all the level crossings of the system in this manner. 
We show in Fig. 2 (b), the morphosis of the eighth wave function around the diabolical 
point (xo,s) = (1/9,0). A noteworthy feature is in the middle right figure (so ~ and xq 
slightly larger than 1/9), where one observes a wave function highly localized between one 
edge and the location of the delta potential. One can generalize this result: By placing a 
strong delta potential at an appropriate location, one can obtain a wave function which is 
localized arbitrary sharply near the edge. 

We next treat the case of two dimension. We consider a quantum particle of mass 
2tt moving inside a rectangle surrounded by the boundaries x — 0, x — R, y = and 
y = 1/R, on which the wave function is assumed to vanish. The eigenstates for this empty 
billiard problem is given by ip nm (x, y) = 2 sin(rmx / R) ■ sm(m7ryR) and their energies e nm = 
7r/4 • {n 2 / R 2 + m 2 R 2 ) (n,m = 1,2...). We place a pointlike scatterer at (x ,y ). Naively, 
we would have the two-dimensional analogue of eqs. (jl]) - However, this is not possible 
since the infinite sum 

g(lu)= jr, ^ m(x ° ,yo) (5) 

n,m=l U ~ £nm 

diverges logarithmically because of the constant density of states per unit energy. Instead, 
we resort to the self-adjoint extension theory of functional analysis. We refer to Ref. 0] for 
a fuller treatment of the problem, and simply start from the results obtained there. The 
eigenvalues of the full problem are given by the solutions of 

G{uj) -s = Q (6) 

where 



DC 



GM = Y. *4«(zo,yo) 



+ 



(7) 



n, m=l 

is the regularized version of G(u>), and s, which we call formal inverse strength, is a real 
number specifying the nature of the pointlike scatterer. The equation (^) can be interpreted 
as the renormalized eigenvalue equation if s is identified as the inverse of the renormalized 



coupling strength. In one dimension, we can consider an analogous equations to eqs.@ 
and(0). In this case, however, each term in the bracket of eq.(|7]) converges separately. This 
allows us to define the inverse strength s in terms of s as a finite quantity. This is the reason 
why we are able to deal with the delta potential in one dimension without any difficulty 
The eigenfunction corresponding to a solution u a of eq.(|5]) is given by 

4>a{x,y)= £ iPnm ^ 0,y °^ (f nm {x,y). (8) 

n,m=l U a ^nm 

We look at the energy of the system as a function of the location of the scatterer, (x , yo) 
with fixed value of s. Because of the mirror symmetry of the system, it is sufficient to 
consider the area xo G [0, R/2] and yo G [0, 1/(2/2)]. In Fig. 3, we show the energy of the 
lowest eleven levels as a function of yo with xo = 0.5R for Fig. 3 (a), and with xo = 0A8R 
for Fig. 3 (b). The value of R is set to 7r/e ~ 1.1557273. The formal inverse strength is 
chosen to s — 1/10. Six diabolical crossings can be spotted in Fig. 3 (including the one 
between the eighth and ninth states which might require further magnification to see the 
avoided crossing off xo = 0.5R). All the crossings involve a state with even n, which is 
unaffected by the presence of the scatterer because of its node line along xo = 0.5R. As a 
result, these states appear flat in Fig. 3 (let us call these flat energies of unaffected states 
e v collectively). 

In general, the level crossing in the system with pointlike scatterer can occur only when 
one of the eigenstates is unaffected by the scatterer. This occurs either along the line 
xq = (M/N)R or along yo = M/(NR) where N and M are the relative primes satisfying 
N > M > 0. There, unperturbed wave functions with (N — 1) nodes in x (or y) direction 
become the solutions of the full problem. The other coordinate yo (or xq) of the crossing 
point is then determined by eq. (||) with uj = e v . Thus we have a systematic way of locating 
all diabolical crossings. As we vary R, the diabolical locations move along xq = (M/N)R 
or yo = Mj (NR). An interesting case is the right square R — 1, where the lowest diabolical 
point (between the second and third states) moves to the center (x ,yo) = (0.5,0.5). We 
observe no sign reversal around this point because this and the other diabolical point at its 
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mirror location y' Q = 1/R — yo merge to form a degenerate diabolical point with the Berry 
phase 2n. 

Going back to the specific case of R = n/e, we display, in Fig. 4, the profile of eigen- 
functions around two of the diabolical points; (a) (xo,yo) = (0-5R, 0.3177/i?) involving the 
second and third states and (b) (x , yo) = (0.5R, 0.3875/i?) for the tenth and eleventh states. 
In case of (a), the third state keeps the unperturbed quantum numbers (n, m) = (2, 1) and 
the energy e v ~ 3.4, while in case of (b), the tenth state has (n, m) = (4,2) and e v ~ 13.6. 
As in one dimension, we learn from these examples that the sign reversal is a result of two 
conflicting ways to connect wave functions with different node structures. One can also look 
at this sign reversal as a rotation of the wave function by 180 degree while the pointlike 
scatterer makes the full turn of 360 degree. As in one dimension, we observe the appearance 
of localized wave functions around the diabolical level crossings. Its possible relevance to 
the "scaring" phenomena [|l(3,[ll]] is an interesting open problem. 



Finally, we view our findings from a perspective and discuss their potential utilities. It 
is important to notice that the rectangular shape of the billiard is a matter of our technical 
choice than the necessary condition for the occurrence of the diabolical crossings. In fact, 
by linearizing G{u) in the neighborhood of e 9 in eq. (|5|), we obtain a second-order algebraic 
equation, with which we can prove the appearance of the Berry phase around diabolical 
points without any reference to unperturbed eigenf unctions. Therefore, we conclude that a 
generic system with a pointlike scatterer possesses diabolical crossings and the associated 
Berry phase. Also, the diabolical crossing discussed here is not a result of the mathematical 
abstraction inherent in a pointlike scatterer. Indeed, we can show through numerical Fourier- 
basis diagonalization that the sign reversal is observed in a billiard with a polygonal obstacle 



of finite size such as the one discussed in Ref. ||12|| . When the size of the obstacle is small, the 



diabolical crossing occurs at a location close to the one predicted for a pointlike scatterer. 
In this case, however, there is no systematic way to locate the diabolical points exactly. It 
is the very pointlike nature of the obstacle that allows the simple expressions such as eqs. 
(§) and (0). We can further hope that the simplicity of the formulation of the pointlike 
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scatterer opens up the possibility of various analytical results concerning the nature of the 
diabolical points such as their geometrical and statistical properties. 

The primary utility of our examples might be a pedagogical one. This is evident when 
we compare them with the earlier examples of the Berry phase in triangular billiards ||13| . 
The rotation in parameter space in our examples literally is the rotation of a pointlike 
scatterer in the coordinate space. Also, the occurrence of the diabolical crossing at the lowest 
possible states including the ground state makes the resulting Berry phase understood in a 
pictorial, intuitive manner. We can also contemplate a possible use of diabolical crossings in 
mesoscopic devices. There may be a situation where an obstacle (or impurity) needs to be 
placed in a certain precise location of the device. This can be achieved by monitoring the 
degenerate resonance energies while moving the obstacle. Another more fanciful use may be 
realized when the phase of wave functions becomes utilized. A quantum state with required 
phase can be obtained by "Berry rotating" an obstacle around its diabolical location. 
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FIGURES 

FIG. 1. Six lowest eigenvalues as function of the inverse strength of the delta potential in one 
dimension: (a) xq = 0.5, (b) xo = 0.47, (c) xq = 0.35, (d) xq = 0.333 and (e) xo = 0.31. 

FIG. 2. The change of the wave function eq. (4) displaying the Berry phase with various 
values of (xq,s) around the diabolic points (a) (xo,s) = (1/2,0) and (b) (xq,s) = (1/9,0) 

FIG. 3. Eleven lowest eigenvalues as function of the perpendicular position of the delta poten- 
tial in two-dimensional billiard while keeping its horizontal position (a) at the center xq = 0.5R, 
and (b) slightly off the center xq = 0A8R. Other parameters are set to R = 1.1557273 and s = 0.1. 

FIG. 4. The change of the wave function displaying the Berry phase around the diabolical 
points (a) (x ,y ) = (0.5R, 0.3177 /R) and (b) (x ,y ) = (0.5R, 0.3875/i?) with R = 1.1557273 and 
s = 0.1. The solid and broken lines indicate positive and negative equilateral lines respectively. 
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